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We present an update to the search for a non-trivial topology of the universe by searching for 
matching circle pairs in the cosmic microwave background [T] using the WMAP 7 year data release. 
We extend the exisiting bounds to encompass a wider range of possible topologies by searching for 
matching circle pairs with opening angles 10° < a < 90° and separation angles 11° < < 180°. 
The extended search reveal two small anomalous regions in the CMB sky. Numerous pairs of well- 
matched circles are found where both circles pass through one or the other of those regions. As 
this is not the signature of any known manifold, but is a likely consequence of contamination in 
those sky regions, we repeat the search excluding circle pairs where both pass through either of 
the two regions. We then find no statistically significant pairs of matched circles, and so no hints 
of a non-trivial topology. The absence of matched circles increases the lower limit on the length 
of the shortest closed null geodesic that self-intersects at our location in the universe (equivalently 
the injectivity radius at our location) to 98.5% of the diameter of the last scattering surface or 
approximately 26Gpc. It extends the limit to any manifolds in which the intersecting arcs of said 
geodesic form an angle greater than 10°. 



I. INTRODUCTION 

The search for a non-trivial topology of the universe 
has enjoyed a long and fascinating history. Using dif- 
ferent methods - from searching for specific topologies 
to the more general circles in the sky approach [2] - 
the cosmic microwave background (CMB) has been an- 
alyzed extensively, looking for any signs that light from 
the same object reaches us by more than one path (see 
[3] for a review of the various suggested methods, includ- 
ing [H |2) I4H22] ). So far, all specialized efforts to detect 
specific topologies as well as the search for matching op- 
posing circles in the sky have failed to detect any sign of 
a non- trivial topology of our universe. 

The circles-in-the-sky method, which we adopt in this 
paper, is based on the following intuitive picture. For 
illustrative purposes, assume that the true topology of 
the universe is a 3-torus, with unit cell size smaller than 
the Hubble horizon (see Figure [I]). This can be thought 
of as a tiling of flat space by identical cubes. An ob- 
server, such as ourselves, performing a series of cosmic 
microwave background (CMB) observations somewhere 
in one of the cubes, has clones identically located in each 
of the other cubes performing the identical series of ob- 
servations. Centered around the observer is the 2-sphere 
of the surface of last scattering at z « 1100, with CMB 
fluctuations imprinted on it. Around the clone of the 
observer on the right, there is another 2-sphere of the 
surface of last scattering. The intersection of both 2- 
spheres is given by a circle. Both observers will look at 
the same ring of temperature fluctuations - albeit from 
different "sides" . Both observers are in fact identical, so 
an observer will see a matching pair of circles: one to 
the left and to the right. Hence, comparing temperature 



fluctuations along circles potentially yields information 
about the topology of our universe. Going away from 
toroidal geometries, it becomes immediately clear that 
the separation angle (the angle between the centers of 
the pair of matching circles) need not be 180°. Also, de- 
pending on the orientability or non-orientability of the 
manifold, circle pairs might have matching temperature 
fluctuations either both going clockwise around the cir- 
cles (non-orientable) or one going clockwise and the other 
anti-clockwise (orientable). 

So far, the search for matching anti-podal circles, i.e. 
circles with separation angles 180°, or nearly anti-podal 
circles [1 , like other topology searches, has only yielded 
lower limits on the size of the Universe, and then only 
for "nearly flat" topologies. 

In this work, we apply the circles-in-the-sky statis- 
tics to searches for circles pairs of all opening angles 
10° < a < 90°, and integer separation angles 11° < 
< 180° with both orientations. This extends the previ- 
ous searches [TJ [6] to cover almost all possible topologies. 
We find what seems to be a systematic effect at two spe- 
cial positions in the sky that produces spurious signals 
for osculating circles. Otherwise, we see no evidence of 
non-trivial topology. 



II. METHODOLOGY 

As already outlined in the introduction, the most in- 
tuitive way to search for a non-trivial topology of the 
universe is by looking for matching pairs of circular tem- 
perature fluctuation patterns in the CMB. 

In order to determine the underlying topology of the 
universe, one would need to scan the full CMB map - 
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FIG. 1: Schematic geometry of the circles-in-the-sky method. 
If the topology of the universe is a torus, the observer, located 
at the center of the 2-sphere of the surface of last scattering, 
sees matching circles of temperature fluctuations on opposite 
sides of the CMB sky. 



at a Healpix 1 resolution of iV s ide = 512 corresponding to 
3 x 10 6 pixels - for matching circle pairs. To conduct 
a search for all topologies, i.e. over all opening angles 
of the circles and all possible separation angles, an enor- 
mous number of circle pairs would need to be analyzed. 
Previously reported searches focussed on nearly flat ge- 
ometries, where matching circle pairs are at almost op- 
posite positions on the sky, i.e. ~ 180° apart (see [HE]). 

Thanks to increasing computing power, we can now 
take a more general approach. In order to obtain an 
acceptable time frame for completing this project, we su- 
perimposed a search grid of N s ^ e = 128 onto the map, 
resulting in 2 x 10 10 circle pairs that need to be compared. 
We compute the goodness of the match for a given circle 
pair on the full resolution, N s ^ e = 512, map. Using a res- 
olution of 512 pixels along a given circle costs about 10 4 
operations per compared pair, leading to a total num- 
ber of operations of 2 x 10 14 per opening angle. On a 
3GHz CPU, this takes about 20 hours for a single open- 
ing angle. Scanning over 200 opening angles corresponds 
to 4000 CPU hours, easily feasible on modern computer 
clusters. 

Thus we search for all possible topologies, albeit on a 
somewhat lower resolution grid. Note that we need to 
run the search twice, once for orientable manifolds and 
once for non-orient able manifolds, i.e. once for circle 
pairs that are oppositely oriented, and once for circle 
pairs that are oriented in the same way. 

A first Ansatz for a circle statistic Sij to measure the 
match between circles z, j would be the convolution of the 
temperature fluctuations along two circles. As discussed 
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FIG. 2: A circle of opening angle a centered at pixel i, and a 
circle centered at pixel i with its phase shifted by /3. 

in [6] , this would lead to a dominance of long wavelength 
(small m) Fourier modes along the circles, making Sij 
rather insensitive to small-scale fluctuations. To com- 
pensate for this, jl] introduced an additional factor of 
the wavenumber m in the convolution (this is equivalent 
to the usual factor of £(£+1) used to scale the two dimen- 
sional power spectrum). Thus, in order to compare two 
circles of opening angle a centered around pixel numbers 
i and j, we employ the circle statistic 
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where f3 is the relative phase between the two circles and 
Ti m (a) is the Fourier transform of the temperature fluc- 
tuation ^ around circle i 

— (i,4>) = T im e— , (2) 

where n = 512 along the circle (see Figure [2|. For a 
search for orientable manifolds, we replace all occurrences 
of T* m (a) in the numerator by T irn (a). Notice that we 
use the conventions of the libfftw package [23 2 : pos- 
itive frequencies are stored in the first half of the ar- 
ray Ti m •> m = • • • § 5 and negative frequencies are stored 
in backwards order from the end of the array, making 
T^ — ^i( n -m) 3 ' We then use as statistic the maximum 
of Sik(a,P) over all relative phases (3 

#™ ax (a) = maxSij(a,P). (3) 



2 |http : / / www . f f tw . org| 

3 Using the real-to-complex and complex-to-real routines enables 
us to save the Fourier transforms for a single circle of a given 
opening angle in an complex array of length ^ instead of n. 
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FIG. 3: Doppler contribution to the matching circle sig- 
nal. Small circles (bottom) have opposite Doppler contribu- 
tions, whereas large circles (top) obtain similar Doppler con- 
tributions. Thus large circles show a smaller deviation from 
S™ ax 1. 



For perfectly matching circles, 5™ ax (a) = 1, and for per- 
fectly uncorrected circles, Sfj a *(a) ~ 0. In practice, 
these ideal values are not realized due to noise contri- 
bution from several different effects. First of all, the 
Doppler effect at the surface of last scattering creates 
a different signal depending on the position of the ob- 
server, making larger matching circles closer to /S™ ax ~ 1 
whereas smaller circles will have 5™ ax < 1 ( see Fig- 
ure [3]). Another contribution comes from line-of-sight ef- 
fects, especially the ISW effect, as the CMB photons from 
different directions traverse through different patches of 
space. The combination of these effects reduces the sig- 
nal in *S^ ax , making potential matching circles less than 
perfect (see Figure 1 in pQ). 

For this analysis, we used the WMAP7 temperature 
maps [24 . Outside the WMAP Kpl2 4 sky mask, we used 
the same template cleaning method as [25]. Inside the 
Kpl2 sky cut, we used the WMAP ILC map 5 . While this 
choice implied that the noise properties and resolution of 
the map differed between the two regions, the effect on 
the circle search was relatively small as the ILC map was 
only used for 5.8% of the pixels. 

In order to estimate the significance of potential spikes 
in the statistics, we created approximate random realiza- 
tions of the CMB 6 . To this end, we compute the a^ m 's 
of the cleaned CMB map. Then, we scramble them 
by randomly interchanging the m index for fixed -Ts, 
and compute a "random" CMB map from the scram- 



4 http:/ / lambda. gsfc. nasa.gov/product/map/drl/intensity_mask.cfm 

5 http://lambda.gsfc.nasa.gov/product/map/dr4/ilc_map_get.cfm 

6 In contrast to the WMAP 1 year data that was examined in 

1 , here the statistics is not noise-dominated. In particular this 
means that the estimator for the false detection rate devised in 
PQ is not applicable for the WMAP7 year data. 



bled a^ m 's. Using this map, we compute the statistics 
for matching (non-) orient able circle pairs for opening an- 
gles a = 20°, 50°, 80° and all integer separation angles 
0. Repeating this 1000 times, we obtain an estimate of 
the probability density function (pdf) of S ma *. While 
there are enough samples to give a reliable estimate of 
the 95.4% (2a) confidence level (CL), the quoted 99.7% 
(3<j) CL is at most a rough estimate. If the sky noise 
were isotropic, then this randomization process would 
generate simulated maps with the same statistical prop- 
erties as the WMAP observations. However, because of 
anisotropics due to spatial variations in the WMAP noise 
(which is larger near the ecliptic plane) and variations in 
the resolution of our map (due to the need to use the 
ILC map in the galactic plane), these simulated maps 
only approximate the WMAP sky maps. 



III. RESULTS 

Implementing the procedure outlined in the previous 
section, we present the results of the searches and de- 
scribe the systematic effects we encountered, both for 
orientable and non-orient able topologies. 

Apart from a possible signal, there are simple random 
statistical fluctuations which are expected to exceed the 
95% CL. To estimate the number of random fluctuations, 
we note that this can be viewed as a series of Bernoulli 
trials with probability p = 0.05 for success, i.e. for a 
spike above 95% CL. Per separation angle we probed 
n — 190 values of the opening angle a (by choosing it 
to lie on a grid deriving from the position of the rings 
in the Healpix scheme for n S id e = 128. The probability 
distribution function for having i excursions above the 
2cr/95.4% CL for a sequence of n Bernoulli trials is given 
by the Binomial distribution 

p l (l-p) n -*(?) , (4) 
whose expectation value is given by 

<iV> = ^ip l (l-p)"-*(j) =np = 8.7. (5) 

Thus, we expect about 9 spikes above the 2a CL per 
separation angle 0. Similarly, we expect about 0.6 spikes 
above the 3a CL per separation angle 0. 



A. Search for Orientable Topologies 

Searching for orientable topologies, we find spikes in 
the distribution of the statistics S as a function of open- 
ing angle a and separation angle (see Figure [4|. The 
red (green) line indicate the 2a(3a) CL as determined 
from scrambling the a^ m 's which we described above. 
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FIG. 4: Searching for orientable topologies: The statistics S (solid blue line) as defined in Equation [I] as a function of opening 
angle a and separation angle (indicated by the numbers in the center of each plot). The red/ green line correspond to the 
2a/3a CL. The spikes come from circles osculating at the galactic anti-center and at / = 109.44°, b — 27.8° (see Figure [7]). 
Removing circles that touch these regions removes all features (see Figure [8}. 
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FIG. 5: Searching for orient able topologies: For separation 
angle = 50°, a plot of the statistics S as a function of 
opening angle a. Without cutting anything, the dashed black 
line shows an extended feature at a = \Q — 25°. Ignoring 
all circle pairs that overlap within 2.5° degrees of the galactic 
anti-center (GAC), the feature mostly disappears (solid blue 
line). The red (green) solid lines show the 95.4% (99.7%) CL 
on S, obtained by scrambling the aim (see the main text). See 
Figure [7] for an explanation of the systematic spikes in S as 
well as Figure [8] for the final result after excising the offending 
circles. 
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FIG. 6: As soon as the separation angle between the two 
circles becomes approximately twice their opening angle a, 
< 2a, the circles start partially overlapping, creating match- 
ing patterns in the statistics S max as for orient able manifolds, 
the "kissing" segments match. 



There is an extended feature whose position (as a func- 
tion of the opening angle a) is correlated with the separa- 
tion angle - it seems to be located at a — \d. Figure [5] 
shows the values of ,S max as a function of a for two circles 
separated by = 50° , with a bump clearly visible around 
a = 25°. This can be understood as the effect of two cir- 
cles partially overlapping (see Figure [6]). If two circles are 
separated by twice their opening angle a they osculate 



("kiss"). For orientable manifolds this means that their 
patterns nearly match up along a segment. This effect 
is almost independent of the absolute magnitude of the 
opening angle a: the fraction of the circles that "kiss" is 
independent of a (up to effects of finite pixel size) . How- 
ever, it is interesting to note that most of the osculating 
circles that give the maximal S for a given separation an- 
gle and opening angle osculate at either of two positions: 
at the galactic anti-center, or at Z = 109.44°, b = 27.8° 
(see Figure [7|). 

Note that there is no known topology which would 
lead to such such a structure of spikes in the S statistics, 
correlating the opening angle with the separation angle, 
confined to just two positions on the sky and appear- 
ing both in the search for orientable topologies as well 
as in the search for non-orientable topologies (see next 
Subsection). Hence it is justified to remove all circles 
that touch either the galactic anti-center or the position 
I = 109.44°, b = 27.8°. The existence of this anomaly 
near the galactic anti-center is perhaps not surprising. 
The anti-center is in the middle of the ILC portion of the 
map, and ILC pixels have much more correlated noise 
properties than the rest of the map. The second point is 
near the region where the noise properties of the maps 
are very non-uniform (see Figure 3 in [26 J. 

First excising circles that osculate at the galactic anti- 
center, Figure^), and then those that osculate at I = 
109.44°, b = 27.8°, Figure [7]i), removes most systematic 
spikes from the S statistics. The peaks that are left for 
separation angles > 170° and opening angle a « 90° 
are caused by the fact that circles with a = 90° have an 
angular diameter of 180°. Thus, two such circles, when 
separated by « 180° will start overlapping, until for 
6 = 180°, they coincide, independent of the topology of 
the universe. 

Hence Figure [8] presents the final result of the search 
for orientable topologies. The number of spikes above 
the 2<j and 3a thresholds are consistent with random 
fluctuations. This leads us to conclude that the data is 
consistent at 99.7% CL with the null hypothesis of no 
non-trivial, orientable topology. 



B. Search for Non- Orientable Topologies 

Performing the search for non-orientable topologies on 
the cleaned CMB map, we again find a multitude of ex- 
cursions above the 3a CL (see Figure [9]). There, we plot 
in blue the statistics S as a function of opening angle 
a and separation angle (the former indicated by the 
numbers in the middle of each panel). The red (green) 
line is the 2a(3a) CL obtained by scrambling the a^ m 's 
as described above. In particular, for low values of the 
separation angle 0, the signal lies consistently above the 
3a CL. The source of these excursions are again circles 
which touch either of two distinct spots: the galactic anti- 
center and I = 109.44°, b = 27.8° (see Figure [l0^). Dis- 
regarding all circles that come within 2.5° of the galactic 
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FIG. 7: For the search for orient able manifolds: Location of the circles pairs with maximal statistics S that lies above the 3a 
CL, colored by opening angle a. Note that the osculating circles at the galactic anti-center and at I = 109.44°, b = 27.8° do 
not hint at a non-trivial topology, a) The highest signal comes from circle pairs osculating at the galactic anti-center and at 
/ = 109.44°, b = 27.8°. b) same as a), but rotated by 180° such that the galactic anti-center is in the middle of the plot, c) 
Removing circle pairs that osculate at the galactic anti-center, the highest signal comes from pairs that osculate at / = 109.44°, 
b = 27.8°. d) Removing circles that osculate either at the galactic anti-center or at / = 109.44°, b = 27.8°, no special position 
on the sky is apparent. 



anti-center, we still find many excursions above the 3<r 
CL (see Figure [Tofr). Finally, also removing circles that 
come within 2.5° of the position I = 109.44°, b = 27.8°, 



there is no structure left (see Figures 10 z and 11). 



Figure [TT] presents the final, negative result of the 
search for a non-trivial, non-orient able topology of the 
universe. The remaining excursions above 2a and 3<r 
are well explained by random statistical fluctuations as 
outlined in the previous subsection. This leads us to con- 
clude that the data is consistent at 99.7% CL with the 
null hypothesis of no non-trivial, non-orientable topology. 



IV. CONCLUSIONS AND OUTLOOK 

We employed the circles-in-the-sky statistics first de- 
vised in [Tj, looking for pairs of matching circles of 
opening angles 10° < a < 90° and separation angles 
11° < < 180°. We positioned the circle centers on a 
grid with N s ^ e = 128, but computed the statistics on the 
full 7V side = 512 CMB map. 

While the WMAP 7 year data brought quite some im- 
provements in the noise of the S^ ax statistics (c.f. Figure 



2 in [T]), we find no hints of a non-trivial topology of the 
universe (see Figures [8] and 11). The new search cov- 
ered a much wider range of possible topologies, and by 
extending the search to circles with opening angles as 
small as 10°, we have extended the previous bound on 
the size of the Universe to 98.5% of the diameter of the 
last scattering surface, or approximately 26Gpc. 

There are systematic effects coming from both mem- 
bers of a circle pair touching either the galactic anti- 
center or the position I = 109.44°, b = 27.8° (see Fig- 
ures [7] and [Toj ) . As these positions appear both when 
looking for orient able and non-orientable manifolds, they 
cannot be of topological origin, but point towards a con- 
tamination of the map at these positions. The galactic 
anti-center region contains significant amounts of galac- 
tic emission. While the ILC maps used in this analysis 
attempt to remove most of this emission, the correlated 
residuals are a likely source of contamination in the circle 
searches. 

We are looking forward to the data release of the 
Planck mission, which will offer an exciting new, sharper 
view of the surface of last scattering, allowing for a better 
search of signs of non-trivial topology by removing noise 
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FIG. 8: Searching for orientable topologies, final result: The statistics S (solid blue line) as defined in Equation [T] as a function 
of opening angle a and separation angle (indicated by the numbers in the center of each plot). The red/ green line correspond 
to the 2a /3a CL. We disregard circles that touch either the galactic anti-center or / = 109.44°, b = 27.8°. The peaks for 
separation angles > 170° and opening angle a « 90° are caused by the fact that circles with a = 90° have an angular 
diameter of 180°. Thus, two such circles, when separated by ~ 180° will start overlapping, until for = 180°, they coincide, 
independent of the topology of the universe. No signs of a non-trivial orientable topology are found. 
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FIG. 9: Searching for non-orientable topologies: The statistics S (solid blue line) as defined in Equation [I] as a function of 
opening angle a and separation angle (indicated by the numbers in the center of each plot). The red/ green line correspond to 
the 2a/3a CL. The spikes come from circles osculating in the galactic anti-center and at / = 109.44°, b = 27.8° (see Figure 10). 
Removing circles that touch these regions removes all features (see Figure 11). 
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FIG. 10: For the search for non-orient able manifolds: Location of the circles pairs with maximal statistics S that lies above the 
3cr CL, colored by opening angle a. Note that the circles that touch either the galactic anti-center or / = 109.44°, b = 27.8° 
do not hint at a non-trivial topology, a) The highest signal comes from circle pairs osculating at the galactic anti-center and 
at I = 109.44°, b = 27.8°. b) same as a), but rotated by 180° such that the galactic anti-center is in the middle of the plot, c) 
Removing circle pairs that touch the galactic anti-center, the highest signal comes from pairs that touch / = 109.44°, b — 27.8°. 
d) Removing circles that touch either the galactic anti-center or I = 109.44°, b = 27.8°, no special position on the sky is 
apparent. 



particularly at smaller separation angles a < 30°. Fur- 
ther advances in computing power will enable a search 
on a full iV S ide > 512 grid of circle positions. 
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